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Abstract 

Upper bounds are obtained for the heat content of an open set D in a 
geodesically complete Riemannian manifold M with Dirichlet boundary 
condition on dD, and non-negative initial condition. We show that these 
upper bounds are close to being sharp if (i) the Dirichlet-Laplace-Beltrami 
operator acting in L 2 (D) satisfies a strong Hardy inequality with weight 
S 2 , (ii) the initial temperature distribution, and the specific heat of D are 
given by 8~ a and respectively, where S is the distance to dD, and 
1< a < 2,1< < 2. 
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1 Introduction 



Let D be a smooth, connected, m- dimensional Riemannian manifold and let 
A be the Laplace-Beltrami operator on D. It is well known (see [11], [14]) that 
the heat equation 

du 

Au=— , x G D, t>0, (1) 
at 

has a unique minimal positive fundamental solution p(x,y;t) where x G D, 
y G D, t > 0. This solution, the Dirichlet heat kernel for D, is symmetric in 
x,y, strictly positive, jointly smooth in x,y G D and t > 0, and it satisfies the 
semigroup property 

p(x,y;s + t) = p(x,z;s)p(z,y;t)dz, (2) 

JD 

for all x,y G D and i, s > 0, where is the Riemannian measure on D. 
Equation (1) with the initial condition 

u(x;0 + ) = ip{x), xeD, (3) 

has a solution 

u^(x;t)= p(x,y;t)ip(y)dy, (4) 
Jd 

for any function ip on D from a variety of function spaces like C\, (D) or LP (D), 
1 < p < oo. Note that u^, G C b (D) if V G C fc (D) or that u^, G L p (D) if 
tp £ L p (D). Initial condition (3) is understood in the sense that (•; t) — > ip (•) 
as t — > + , where the convergence is appropriate for the function space of initial 
conditions. For example, if ip G Cb (D) then the convergence is locally uniform, 
or if ip G i p (Z?), 1 < p < oo then the convergence is in the norm of LP (D). In 
general, (4) is not the unique solution of (l)-( 3). However, it has the following 
distinguished property: if ip > then is the minimal non-negative solution 
of that problem (and if ip is signed then = — u^,_). If D is an open 
subset of another Riemannian manifold M and if the boundary dD of D in M 
is smooth then the minimality property of u^, implies that, for any t > 0, 

lim u^,(x; t) = 0. (5) 

x—tdD 

If dD is non-smooth then (5) can still be understood in a weak sense. Expression 
(4) makes sense for any non-negative measurable function ip on D, provided the 
value +oo is allowed for u^. It is known that if G L\ oc (D x M + ) then is 
a smooth function in D x R + and it solves (1) (see p. 201 in [14]). For any two 
non-negative measurable functions ipi,ip2 on D, we define for t > 

Qin^it) = // p(x,y,t)ip 1 (x)ip 2 (y)dxdy. (6) 

J JDxD 

Using the properties of the Dirichlet heat kernel we have for < s < t 

Qi>ui>2(t) = / u^(x; s)u^, 2 {x; t - s)dx. (7) 

Assuming that D is an open subset of a complete Riemannian manifold M, 
Qip 1 ,ip 2 (t) nas the following physical interpretation: it is the amount of heat in 
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D at time £ if D has initial temperature distribution ipi , and a specific heat ip2 , 
while the dD is kept at fixed temperature 0. 

This function has been subject of a thorough investigation. Its asymptotic 
behavior for small £ is well understood if D has compact closure with C°° 
boundary, and both tpi and ip2 are C°° on the closure D of D. In that case 
Qi/>i ,i/> 2 (£) has an asymptotic series in £ 1//2 , and its coefficients are computable 
in terms of local geometric invariants [2, 12]. No such series are known if D is 
unbounded, or if either the initial data or dD are non-smooth. 

In this paper we will obtain upper bounds for the heat content Q^^ip) 
under quite general assumptions on D and on ipi and ip2- 

We are particularly interested in the situation where D is a open subset of 
another manifold M, and where tpi (x) and ip2 (x) blow up as x — > dD. In order 
to guarantee finite heat content for £ > 0, sufficient cooling at 3D needs to take 
place. This will be guaranteed by a condition on D, that is formulated in terms 
of a Hardy inequality. Note that in this setting Q^i ,V>2 (*) ma y be unbounded 
as £ — > + , and one of the interesting points of this study is to obtain the rate 
of convergence of <3i/>i,V>2 (£) to +oo as £ — > + . 

Given a positive measurable function h on a manifold D, we say that the 
Dirichlet Laplacian acting in L 2 (D) satisfies a strong Hardy inequality with 
weight h if, for all w £ C c °°(-D), 

Here, and in what follows, we put f D f = j D f(x)dx if this does not cause 

confusion. We also put \D\ = j D 1, and ||/|| p = ( J D \f\ p ) 1 ^ P - A typical example 
of a Hardy inequality is when D is an open subset of another manifold M, and 

h(x) = c 2 S(x) 2 , (9) 

where c > 2 is a constant, S is the distance to the boundary, 

S(x) = mm{d(x,y) : y £ dD}, 

and d(x,y) is the geodesic distance from x to y on M. Both the validity and 
applications of Hardy inequalities with weight (9) have been investigated exten- 
sively [1], [7], [9], [10], [11], [4]. For example, inequality (8) holds with weight 
(9) with c = 4 if D is simply connected with non-empty boundary in ffi 2 , with 
c = 2 if D is convex in R m , and for some c > 2 if D is bounded with smooth 
boundary in R m . 

In [3] it was shown that if D has finite measure and satisfies the Hardy 
inequality with weight h, and if ip is a non-negative measurable function on D, 
such that, for some q > 1, 

<oo, (10) 

then, for all £ > 0, 

Q*,i(t) < (iT^T)) 17 " - u, {■;t)\\\ /q t- 1/ \ (11) 
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where Q\ t \ is defined by (6) for tpi = ip 2 = 1, that is, 

Qi,i(*) = / Mi(ar;t)da;= // p(x,y;t)dxdy. 
jd JJdxd 

A similar estimate holds for arbitrary open sets D C R m , satisfying the 
Hardy inequality with weight h. If t\) is a non-negative measurable function on 
D such that, for some q > 1, 

||max{V,l}/i 1/9 || 9 / (9 -i) <oo, (12) 

then, for alH > 0, 

Q*,i(t) < aWII^ 1 ^!!,/^-!)!^ 1 /^*- 1 "!!,*- 1 ^ 9 - 11 , (is) 

where; 

/ „ \ (2?-l)/(9(9-l)) 
- A" 1 /? ( _J_ 



«(<?)= 4 • ( 14 ) 

Below we give a sufficient condition for the fmiteness of Qvi.V^M f° r au 
i > 0, and reduce the problem of finding upper bounds for Q^.^t) to the case 

"01 = ^2- 

Theorem 1. Le£ ipi and ip2 be non-negative and Borel measurable on a manifold 
D. 

(i) If Qii>i,4>i{t) < oo, i — 1,2, for all t > 0, then Q^ 1 ^ 2 {t) < oo for all t > 0, 
and 

Q^M*) ^ (*)0^i,*2(*)) 1/2 . f > °- ( 15 ) 

(ii) If Q$ it i(i) < oo, i — 1, 2, /or aZZ i > 7 and i/ 

c t := swp p(x,x;t) < oo, t > 0, (16) 

x£D 

then 

Qtufoit) < c t / 3 Q^ u i{t/3)Q 4 , 2A (t/3) < oo, t > 0. 

Our main results are the following three theorems, in which we assume that 
D is a Riemannian manifold that satisfies the Hardy inequality with some weight 
h, and ip is a non-negative measurable function on D. In particular we do not 
assume any smoothness conditions on dD, nor do we assume that D has finite 
measure or that D is bounded. 

Theorem 2. If \D\ < oo, and if there exists 1 < q < 2 such that 

< oo, (17) 

then, for allt>0, 

(4-<j)/<j 

Q*At) < (2(g - i))2/ g H^II^-dII 1 - »i (-:*)lli 2 " g)/9 t- 2/g - (18) 
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Theorem 3. If 1 < q < 2 is such that (17) holds and that 

\\h 1/q \\ q / {q -i) < oo, 

then 

Q*At) < HQ) \\Tph 1/q \\ 2 q / {q -i) ll^ll^)'" 1 ^- 1 ^'-^, * > 0, (19) 
w/iere 

6(g) = 2( 4 - 3 «)/(«(«- 1 » f^-yj . (20) 

Theorem 4. If0<r<2, and 1 < q < 2 are such that 

W\\g < 

and 

ll^-^llfa-D/^oo, 

Qv^W < ll^ll,!^ 2 -^ 17 '!!,/^-!)*- 1 ^, t > o. (21) 

In Theorem 5 in Section 3 we use the bounds of Theorems 2 and 4 together 
with (15) to obtain an upper bound for the heat content of D, when D satisfies 
a Hardy inequality with weight (9), and ipi(x) = S(x)~ a and ip2{x) = $(x)~^, 
where a, (3 £ (1, 2). Even though the bounds in e.g. 2 and 4 look very different, 
both of them are needed to cover the maximal range of a and /? in Theorem 5. 

Theorem 2 has a curious consequence. We claim that if a manifold D has 
finite measure \D\, and is stochastically complete then no Hardy inequality holds 
on D (which confirms the philosophy that the Hardy inequality corresponds to 
cooling that comes from the boundary). Indeed, stochastic completeness means 
that mi = 1. In this case, ||1 — ui(-;t)|| 1 = so that we obtain from (18) that 
Qip,i{> (t) — whenever function ip satisfies the condition (17) for some q e (1, 2). 
However, if h is finite then it is easy to construct a non-trivial function tp that 
satisfies (17): choose any measurable set S with finite positive measure such 
that h is bounded on 5, and let ip = Is- Then (17) holds with any q > 1 while 
Qip,4> {t) > so that we obtain contradiction. Of course, without the finiteness 
of \D\. the Hardy inequality may hold on stochastically complete manifolds like 
R m \ {0}. 

This paper is organized as follows. In Section 2 we will prove Theorems 1, 
2, 3 and 4. In Section 3 we will state and prove Theorem 5. Finally in Section 
4 we obtain very refined asymptotics in the special case of the ball in R 3 with 
■01 (x) = 5(x)~ a ,a < 2,tp 2 (x) = 5(x)-f\f3 < 2, and a + /3 > 3 (Theorem 7). 
This special case shows that the bound obtained in Theorem 5 is close to being 
sharp. Moreover it suggests formulae for the first few terms in the asymptotic 
series of a compact Riemannian manifold D with the singular data above. 
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2 Proofs of Theorems 1, 2, 3 and 4 



Proof of Theorem 1. In both parts, it suffices to prove the claims for non- 
negative functions ip2 from L? (D). Arbitrary non- negative measurable func- 
tions -01 * "02 can be approximated by monotone increasing sequences of non- 
negative functions from L 2 (D), whence the both claims follow by the monotone 
convergence theorem. 

To prove part (i) we use symmetry and the semigroup property, and obtain 
by (7) for s = t/2 that 

<9v>i,V>2(*) = / u^, 1 (x;t/2)u^ 2 (x;t/2)dx 
J D 

< «& (x; t/2)<tej u% 2 (x; t/2)<bj 

1 /2 

= (Qii>i,ii>i(t)Qii>2,fo(i)) 

It follows from (2) and (16) that 

p(x, y; t) < {p{x, x- t)p(y, y; t)) 1 ' 2 < c t . (22) 
To prove part (ii) we have by (22) that 

p(x,y;t) = // p(x,z 1 ;t/3)p(zi,Z2;t/3)p(z2,y;t/3)dz 1 dz 2 

J J DxD 

< c t /3 M i( a; ;i/3)Mi(2/;t/3). (23) 
This together with definition (6) completes the proof. □ 

For the proofs of Theorems 2, 3, 4, choose a sequence {D n } that consists of 
precompact open subsets of D with smooth boundaries such that D n C D n+ \ 
and \J n D n = D. Obviously, Hardy inequality (8) remains true in any D n with 
the same weight h, because (D n ) C (D) . Moreover, we claim that 
(8) holds for any function w E C (D n ) n C 1 (D n ) that satisfies the boundary 
condition w|a_o„ = 0. Indeed, if J D \\7w\ 2 = oo then (8) is trivially satisfied. 
If J D |Vw| 2 < oo then w belongs to the Sobolev space W 1 - 2 (D n ). Extend 
function w to D n+ \ by setting w = in D n+1 \ D n . Due to the boundary 
condition w|an„ = 0, we obtain that w n £ W 1 ' 2 {D n+1 ). Since w is compactly 
supported in D n+ i, it follows that w £ W 1,2 (D n+ i) where W 1,2 (ft) is the 
closure (f2) in W 1 ' 2 (Q). Since the Hardy inequality (8) holds for functions 
from C£°(D n+ i), passing to the limit in W 1 ' 2 {D n+ {) and using Fatou's lemma, 
we obtain that w also satisfies (8). 

Assume for a moment that the statements of the theorems have been proved 
in each domain D n . Then one can take the limit in (18), (19), (21) as n — > oo, 
and obtain the statements for D. Indeed, the left hand side of these inequalities 
is Qfy (t) = ff D ^ xDn p Dn (x,y;t)ip(x)'ip(y)dxdy, wherep D „ is the Dirichlet heat 
kernel for D n . This converges to ^ (t) as n — > oo. The right hand sides of 
(18), (19), (21) contain various L p (D n )-norms that can be estimated from above 
by the LP (D)-norms. The only exception is the term ||1 — J D po n (■, y, t)dy\\\ 
in (18) that is decreasing as n — > oo. If \D\ < oo then 1 £ L 1 (D) so that the 
passage to the limit is justified by the dominated convergence theorem. 
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Hence, it suffices to prove each of the statements for D n instead of D. Re- 
naming D n back to D, we assume in all three proofs that D is a precompact 
open domain with smooth boundary in M. 

Another observation is that all inequalities (18), (19), (21) survive the in- 
creasing monotone limits in if). So it suffices to prove them when is bounded 
and has a compact support in D, which will be assumed below. Furthermore, 
since all the statements of Theorems 2, 3, 4 are homogeneous with respect to 
ip, we can assume that < tp < 1. If ip = then there is nothing to prove; 
hence, we assume that ip is non-trivial. Then (x; t) is smooth and bounded 
in D x (0, +oo) and positive in D x (0, +oo). 

Proof of Theorem 2. Let v be the outwards normal vector field on dD. Using 
the Green's formula, we obtain 

d f q f q-ldUijj 

- ' u l = -i / «3 - - 



dtJ D u + *J»* dt 

Jd 



dD v dv j D 



(Vuf\V^) 

= q(q-l) f <" 2 |V^| 2 , (24) 
J D 

where we have used that q > 1 and, hence u 1 ^ 1 = on dD. Observing that 



■ q/2 eC(D)nC 1 (£>), 



2 = ^„r 9 ivu-' a 



4 -v i vu v-i » 

and applying the Hardy inequality (8) to u q / 2 , we obtain that 
By Holder's inequality we have that 



q \ 1/9 



f ll^l^-i)- (26) 



By (25) and (26) we conclude that 

-i["l> ^^^ hl/ %/( q -i) ■ (27) 



dtJ D *~ q 

Note that the function t i->- Q^.-^(t) = \\u^ (■;t/2)\\ 2 is decreasing in i, which, 
for example, follows from (24) with g = 2. Integrating differential inequality 
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(27) with respect to t over the interval [t, 2t] gives that 

J D n% > fcllll^H-^^ (Q^(2t))H. (28) 
On the other hand, using 1 < q < 2 and the Holder inequality, we obtain 

that is, 

/ < (*)) 2 "' (Q^(2*)) a-1 ■ (29) 

Combining (28) and (29) yields 

QtAW < 4^3i)ll^ 1/a H2/(,-i) WwW) 2 "' 4 " 1 - ( 3 °) 

Estimating Q$,i by (11), we obtain 

Q*,m < ( 4(g _ i))2/ g II^I^-dII 1 ~ "i(-;*)lli 2 - g)/g *- 2/g , 

which completes the proof. □ 

Proof of Theorem 3. Since ^ < 1 we have that (12) is satisfied. We obtain by 
(13) and (30) that 

This completes the proof of Theorem 3 since, by (14) and (20), 
2l/( ' _1) 4(^31) °(*) a "' = ft («)- 

□ 

Proof of Theorem 4- By the arithmetic-geometric mean inequality, we have 

V>(*)V>(y) < ^ ^{x) r ^{y) 2 - r + ^{x) 2 -^{ y y) . 

By non-negativity and symmetry of the Dirichlet heat kernel 

Q^rt (t) < ! u^ 2 - r . (31) 
Jd 

Next, Holder's inequality yields 

j D u ^ 2 ~ r < (j D u l'-l) 1/q W^^h/^y (3 2 ) 



By (25) we have 



A [ u * > /" u * I (3 3) 

dtJ D u r- q J D u ^h- (33) 



Combining (31), (32), (33) we obtain that 

«♦*<*»* £ -loFTj! (/>) i* 2 -*"'^.)- 

Since the function t i— > Q ,/,,,/, (i) is decreasing in t, we obtain by integrating the 
differential inequality (33) with respect to t over the interval [0, t] that 

and (21) follows. □ 

3 Singular initial temperature and singular spe- 
cific heat 

Below we make some further hypothesis on the geometry of D, and obtain an 
upper bound for the heat content for a wide class of geometries using Theorems 
2 and 4, and (15), if the initial temperature distribution and specific heat are 
given by <5~ Q , 1 < a < 2, and S^ 13 , 1 < (3 < 2 respectively. 

Theorem 5. Let D be an open set in a smooth complete m- dimensional Rie- 
mannian manifold M , and suppose that 

i. The Ricci curvature on M is non-negative. 

ii. For x G D, 

1>a(x) = S(x)~ a . 

Hi. D has finite inradius i.e. pn — sup{S(x) : x G D} < oo. 

iv. There exist constants kd < oo, d G [m — l,m) such that 

[ 1 < n D p m - d , 0<p<p D . (34) 

J{xeD:S(x)<p} 

v. The strong Hardy inequality (8) holds with (9) for some c > 2. 
If 1< a < 2, 1 < /3 < 2, andife>0 then 

(*) = 0(t- £ +( m - d -«-«/ 2 ), t -> 0. (35) 

Proof. Note that (hi) and (iv) in Theorem 5 imply that \D\ < K,up r ^, d < oo. 
By (15) it suffices to prove (35) in the special case a = fi with 1 < a < 2. In 
order to estimate ||1 — ui(-;i)||i in Theorem 2 we rely on the following lower 
bound for u\ (Lemma 5 in [5]). 

Lemma 6. Let M be a smooth, geodesically complete Riemannian manifold with 
non-negative Ricci curvature, and let D be an open subset of M with boundary 
dD. Then for x G D,t > 

Ml (x;i)>l-2( 2+m )/ 2 e -^ 2 /( 8t ). 
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To prove (35) we first consider the case 

(2 + m- d)/2 < a < 2. (36) 
This set of a's is non-empty since d G [m — 1, m). By (9) we have that 

lllM 1/9 H,/(,-i) = c 2/9 [JjV-vVto- 1 )) . (37) 

Denote the left hand side of (34) by ujd(p)- Then we can write the right hand 
side of (37) as 

c 2 ^U pV-<rilto-Vw D {dp)\ . (38) 

An integration by parts, using (36) shows that (38) is finite for 

2 - m + d . . 

q < — . 39 

a — m + d 

Since a satisfies (36), we have that the right hand side of (39) is in (1,2). We 
now choose e > such that 

2 - m + d A i + ')- 1 e(l,2), (40) 

a — m + d 

and choose q equal to the left hand side of (40). By Lemma 6 and (34) we have 
that for t ->■ 



||l-ui(-;i)||i = / {I - u(x;t))dx 

J D 

K 2 (m+2)/2 f e -S 2 /(8t) 
JD 

< 2 (m+2 ^ 2 f e~P 2 ^uj D (dp) 

= 2 (™+2)/2 e -pL/(8t)| £ )| + 2 (m-2)/2 Kz?t -l f'" p m-d+l ^/(St) df) 

Jo 

= 0{& n - d ^ 2 ). (41) 

By Theorem 2 and (37)-(41) we find that for all a satisfying (36) and all e > 
satisfying (40) 

Q^.V-a (*) = 0{t-< a - m+d ^ rn - d - 2a ^ 2 ), t -> 0. (42) 

We conclude that (35) holds for all a = satisfying (36), and all e > 0. 
Next consider the case 

1< a < (2 + m-rf)/2. (43) 
This set of a's is again non-empty since d E [m — 1, m). By (34) we have that 

W\\ q = [J LJ D (dp)p- ar ^ ^ < oo (44) 
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for 

arq < m — d, (45) 

and 

ll^ 2 -^ 1/9 ll 9 /( 9 -i) = (^ + WD(dp)pl 2 - a < 2 - r Wl«-»y ^ < oo (46) 

for 

«g(2-r)-2 <m _ d 

q- i 

The optimal choice for r is henceforth given by 

r = 2{aq-l)a- 1 q- 2 . (48) 

By (43) we also have that a > 1. Hence r e (0, 2). The requirements under (45) 
and (47) become with this choice of r that 

q < 2{2a + d-m)- 1 . (49) 

Since a satisfies (43), the right hand side of (49) is in (1,2). We now choose 
e > such that 

2((2a + rf-m)(l + 2e))- 1 e (1,2), (50) 

and choose q equal to the left hand side of (50). By Theorem 4 and (44)- (49) 
we find that for all a satisfying (43), and all e > satisfying (50) 

Q^ a (i) = (t-< 2a - m+d ^ m - d - 2a ^ 2 ), t -> 0. (51) 

We conclude that (35) holds for all a = f3 satisfying (43), and all e > 0. 

To prove (35) for the limiting case a — (3 = (2 + m — d)/2 := a c we note 
that Q$,<t>{t) is monotone on the positive cone of non- negative and measurable 
ip and <j). Let a = a c + e where e is such that a £ (a c , 2). Since 

< P a D a ^ a - 

we have by (42) that 

_ Q^-e(2+e+(d-m)/2)+(m-d-2a c )/2^ ^) 

We conclude that (35) holds for a = (3 = a c , and all e > 0. □ 

4 The special case calculation for a ball in R 3 

In this section we show by means of an example that the upper bound obtained 
in Theorem 5 is close to being sharp for a < 2,(3 < 2, a + (3 > 3. 
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Theorem 7. Let B a = {x e M 3 : \x\ < a}. If a < 2, [3 < 2, a + [3 > 3, J 6 N 
f/iere exist coefficients bo, &i, • • • depending on a and j3 only such that for 
t -> 

<3^,^ (*) = 47rc a ^a 2 i( 1 - a -' 3 )/ 2 - 47r( Ca _ li ^ + Ca^at*- 2 -"-™ 2 

+ 47rc a _ 1;/3 _ 1 t( 3 - a -' 3 )/ 2 + ]T b ja 3 -*- a -W 2 + 0(t( J+1 )/ 2 ), (53) 

j=0 



where 



c a3 =2-«-^- 1 /2r((2- a -/3)/2) 

x /V" + ^((l - p)^" 2 - (1 + p) a+ ^ 2 )dp, 
Jo 



(54) 



and 



6 = -87r((a + (3- l)(a + /3 - 2)(a + /3 - 3))"\ 

61 =0, 

6 2 = 87ra/3((a + /3 + l)(a + /3)(a + (3 - 1))"\ 

6 3 = 0. (55) 

We see that the leading term in (53) jibes with (35) since (9) holds for some 
c > 2, and (34) holds with d = m — 1. 

Theorem 7 suggests that for any precompact D with smooth dD in M, and 
for a < 2, /3 < 2, a + (3 > 3 and t -> 



t (2-a-P)/2 



Q^eit) =c a ,p ( t {1 - a -' ))/2 -2- 1 {c a - 1 .f j +c a j 3 - l ) [ L gg 

JdD JdD 

+ [ (c 1 L gg L hh + c 2 L gh L gh )t^- a -^' 2 + 0(1), (56) 
where ci and C2 are constants depending on a and (3 only, and which satisfy 

4ci + 2c 2 = Ca-1,/3-1, 

and where L gs is the trace of the second fundamental form on the boundary of 
3D oriented by an inward unit vector field. Since J gB 1 = 47ra 2 , J gB L gg = 8wa 
and J gB (ciLggLhh + C2L g hL g h) = 16nci + 871x2, we see that (56) holds for the 
ball in K 3 . 

The proof of Theorem 7 rests on the following result (pp.237, 367-368 in [8]). 

Lemma 8. Let B a as in Theorem 7, and let the initial datum be radially sym- 
metric i.e. ipi(x) = f(r), where r = \x\. Then the solution of (1), (3), (5) is 
given by 



X-t) = (Mr 2 )- 1 / 2 /%'/(/) ]T(e- (2rla - ^TA^) _ e -(2na+r+r'f/m )dr ^ 
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To prove Theorem 7 we have by Lemma 8 that 
Q^M ^{An/tf 2 f f rr > {a _ r) -a (a _ r /)-/3 



X ^( e -(2™-r+r') 2 /(4t) _ g _ (2„a+r+r') 2 /(4 t ) ) (57) 

where 5 a = [0, a] x [0, a]. Substitution of a — r = p and a — r' = q in (57) gives 
that 

Qi> a ^fi (*) = + Ai + A 2 + B, 

where 

A = {^/t) l ' 2 a 2 J I p ->-^( e -(P-9) 2 /(4t) _ e -(P+^ 2 /m)dpdq, 
A l = -{^/t) l ' 2 a JJ (p + g)p- a g-^(e-^-«) a /W - e-<*-«> a /< 4t >)dpdg, 

A 2 = (4^A) V2 y y p 1 -«g 1 -/3( e -(f-9) 2 /(4t) _ e -(jH- fl ) a /(4*)) dpdgj 

and 

B = {Air/t) 1 ' 2 // (a - P)(a - q)P~ a q~ J2(e- {2na+p - q)2/{it) 

+ e -(2na+q-p) 2 /(At) _ e -(2na+q+pf / (it) _ g - (2na-g-p) 2 /(4t) 

(58) 

We have the following. 

Lemma 9. //l<a<2,l</3<2 f/ien for t ^ 

B = -Stt 1 /^- 1 ^-' 3 * 3 / 2 + 0(t 2 ). (59) 
Proof. The integrand in (58) can be rewritten as 
(a - p)(a - q)p~ a q-P ^ e -(2«a- P - 9 ) 2 /(4 t ) 

n>l 

x ^ e (p- 2 " a )«/* _|_ e (9-2n<»)p/t^l _ g-P9/*) _ ^ _ g-2pna/t^ _ g -2gno/t^ 

(60) 

The contribution from the terms with n > 2 in (60) is bounded in absolute value 

by 

2ay-V" /3 i- 1 E e- a2( "- 1)2/ *(l + 2nVt- 1 ). 

n>2 

After integrating with respect to p and 9 we see that this term contributes at 
most 0(e- a2 /( 2t )) to B. Next we will show that the main contribution from the 
term with n = 1 in (60) comes from a neighbourhood of the point (p, q) — (a, a). 
Let 

Ci(a) = {(p,q) eM 2 : a/3 < p < a,a/3 < q < a}, 
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and 

C 2 (a) = S a \C 1 (a). 

On C2(a) we have that 2a — p — q > 2a/3. Hence the term with n = 1 in (60) 
is bounded on 02(a) in absolute value by 

2(o - p)(a - g)? 1 -"? 1 ^-^-" 2 /' 9 ^! + 2a 2 r 1 ). (61) 

Integrating (61) over 6*2(0) gives a contribution which is bounded by 
0(e~ a /( 18 *)). In order to calculate the contribution from the term with n = 1 
on C\ (a) we use the expression under (58) instead. First we note that 2a+p—q > 
2a/3, 2a + q — p > 2a/3, 2a + p + q > 8a/3. Hence the first three terms in the 
summand of (58) with n = 1 give after integration over C\(a) a contribution 
0( e - a2 /(i8*)). Putting all this together gives that 



B = - (4 7 r/t) 1 /2 ff ( a _ p)(a _ q ) p -° q -P 

J J Ci(a) 

X e-^ a -i-rfl^)dpdq + 0( e - a2 /< 18t )). 

Noting that 

p- a q- f> = a- a - p + 0(a - p) + 0(a - q) (62) 
uniformly in p and q yields after a change of variables that 

B = -(^/t) l ' 2 a- a -' 3 ff 

JJS a/ 

X (1 + 0(p) + 0(q))dpdq + 0(e- a2 /( 18 *)), 

which agrees with the right hand side of (59). □ 

By taking higher order terms of the form (a — p) ni (a — q) n2 in (62) into ac- 
count one can determine the coefficient t^ +3 ^ 2 ,j = 0, 1, 2, • ■ • in the expansion 
of B. 

To complete the proof of Theorem 7 we rewrite A , A x and A 2 respectively 
as follows. 

A =(47r/t) 1/2 a 2 / dp dq + dq dp) 
\Jo Jo Jo Jo J 

x p- a q-l 3 (e- ( - p - q)2/{4t) - e -(P+9) 2 /(4*)) 

=(4 7 r/t) 1 / 2 a 2 f p^-Pdp f (p- a +p-P) 
Jo Jo 

x(e -p 2 (l-p) 2 /(4*)„ e ^ 2 (l+p)V(4 t ) )dp 

=4 7 ra 2 c Q ,^ 1 - a -«/ 2 

- (47r/i) 1/2 a 2 ^ p^-f'dp f(p- a +p- tj ) 

J a Jo 

x (e -f 2 (i-P) 2 /(4*) _ e -p 2 (i+ P r/m )d ^ (63) 
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POO 

A x =- 47ra(c a _ li/3 + c^^ 2 ^-^' 2 + (Air/tf^a / p 2 - a - f3 dp 

J a 

x f dip 1 -" + P - a + p 1 -? + p -^)( e -p 2 (i-P) 2 /(4 t ) _ e -p 2 a+p) 2 /m )dp ^ 

Jo 

(64) 

and 

(-OC 

^2 = 4^C ct _ 1 ,^ 1 ^ 3 - Q -«/ 2 - (47TA) 1 / 2 / P 3 -«-^p 

•/ a 

x f 1 dip 1 -* + p i-^)( e -p 2 (i-p) 2 /(4 t ) _ e -p 2 a+p) 2 /m )dp . (65) 

Jo 

The terms to be evaluated in (63), (64) and (65) are all of the form 

(4^A) 1/ V~-' / p 1+ 3- a - f3 dp p-T( e -P 2 (l-p) 2 /(4t) _ e -P 2 (l+P) 2 /W )d/3j 

(66) 

where j = 0,1,2 respectively. Following arguments similar to the proof of 
Lemma 9 we see that the contribution of the integral with respect to p e [0, 1/2) 
in (66) is at most <3(e~ a /( 18 *)). Furthermore 

fCO fl 

(4n/t)^ 2 a 2 ^ p 1+ i- a -f 3 dp p-<e-r 2(1 +rf/^dp = 0(e- a2 /^). (67) 

Ja Jl/2 

Hence the expression under (66) equals 

fOO fl 

{^/ty/ 2 a 2 -n p 1+ i- a -Pdp p-^e-^^-rf/^dp + Oie-^l^). (68) 

Ja Jl/2 

Expanding p^ 1 about p = lwe obtain that 

\p-i - 1 - 7 (1 - p) ~ 2^7(7 + - P? 

- 6-^(7 + 1)(7 + 2)(1 - P f\ < C{1 - p)\ < p < 1/2, (69) 

where C depends on 7 only. By (69) and (68) we obtain that (66) is equal to 

2^(a + l)- 1 ^-*-' 3 + 47T 1 / 2 7 ( a + - j)- l a 2 - a -' 3 t 1 ' 2 

+ 27r 7 (7 + l)(a + P-J + l)' 1 ^-"-^ 

+ 87r 1 / 2 3- 1 7 ( 7 + 1)( 7 + 2)(a + (3 - 3 + 2)- 1 a- a ~H s ' 2 + 0(t 2 ). (70) 

It remains to compute the coefficients b$,bi an d bi in Theorem 7. Altogether 
there are eight terms which contribute to the terms in (70): 

3=0, J = a, 7 = P 

3 = 1, 7 = a - 1, 7 = £ - 1, 7 = 7 = /? 
j = 2, 7 = a - 1, 7 = /3 - 1 . 

Summing these eight terms yield the expressions for bo,bi and 62 under (55). 
To calculate 63 we have that the above eight 7(7 + 1) (7 + 2) terms in (70) cancel 
the contribution from (59). This completes the proof of Theorem 7. 
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